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1. Introduction 

The subject of this paper is a class of two- loop electroweak contributions 
to the muon g — 2 containing a fermion triangle along with a virtual photon 
and Z boson, as shown in Fig. 1. In the standard model all charged fermions 
contribute to the triangle loop. Individual contributions of fermions lighter 
than the Z boson are enhanced by large logarithms ln(Mz/m), where m 
denotes the mass of the fermion in the loop or of the muon, whichever is 

* Based on talks given by AC or AV at the XLIII Cracow School of Theoretical 
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clear Physics (New York, May 2003), and the International Symposium on Lepton 
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Fig. 1. Effective Z-f-f* coupling induced by a fermion triangle, contributing to 



EW 



heavier. Those large logarithms were first found in [1] where the diagrams 
in Fig. 1 where evaluated only for leptons in the triangle loop. However, 
it was pointed out in [2] and explicitly shown in [3] that such logarithms 
cancel in sums over all fermions of a given generation (as long as m < Mz 
for all fermions in the generation). The source of this cancellation can be 
traced back to the cancellation of anomalies, given by the same fermionic 
triangles, within a given generation. It is required for renormalizability. 

Some subtlety here is that the anomaly cancellation refers to the lon- 
gitudinal part of the axial Z boson current, while both transversal and 
longitudinal parts of the fermionic triangle contribute to g — 2. At the level 
of free quarks the cancellation is, of course, explicit, so the only question is 
whether it can be spoiled by strong interactions. Since In Mz arises from the 
loop momenta much larger than the hadronic scale Aqcd perturbative QCD 
should provide the answer. For the longitudinal part the Adler-Bardeen the- 
orem [4] guarantees the absence of gluonic corrections governed by a s {Mz)- 
We will describe below why this is also true for the transversal part con- 
firming the cancellation of In Mz within a generation. 

For the light quark loops strong interactions become essential at the 
range of loop momenta on the order of hadronic scale. Not only is a s ~ 1 
but also nonperturbative phenomena defining hadron masses are crucial in 
this range. Namely, these dynamical phenomena — not the current masses 
of light quarks — define effective infrared cutoff in the logarithmic integrals. 
While an exact calculation in this range is not possible, in the paper [3] 
(CKM) a crude approximation was used: "constituent masses" of 300 MeV 
assigned to quarks u, d and 500 MeV to the quark s played the role of the 
infrared cutoff. The contribution of the first two generations of fermions 
was found to be 



R . (m u m c ) 4 /3 49 8tt 2 
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~ -8.3 • 10" 11 (with m c = 1.3 GeV). (1) 
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Some uncertainty was assigned to this result, but the cancellation of In Mz, 
present in that model, was believed to be rigorously valid. 

Around the same time the paper [5] (PPR) appeared, in which the 
fermion triangle contributions of Fig. 1 were also studied. The result of 
the PPR paper for the first two generations was 



,PPR/„ „. a. ,. „ „\ . . Gf m l a 



Aa^ (e,u,d;v,c,s) = ^g^~ 

~ -8.7 -lO -11 . (2) 
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This calculation was done in the chiral limit, so that the u,d,s masses were 
replaced by in (2), just like the electron mass. The Mz dependence of 
the logarithms would cancel if the first coefficient (—14/3) were replaced 
by —12/3 = —4. We have argued that —4 is indeed the correct coefficient 
and that the factor —14/3 arises due to an incomplete accounting of u, d, s 
contributions in Eq. (20) in [5]. 

Since the numerical values of (1) and (2) are very similar, and since 
analytical results in [5] were given for individual flavors and not for their 
sum, the presence of a residual In Mz in the final result of [5] was incon- 
spicuous. The difference was first pointed out by Mingxing Luo in a private 
communication to W. Marciano. 

More recently, the authors of [5] together with M. Knecht have revisited 
that issue in a detailed study [6]. They maintained the finding of [5] and 
argued that the cancellation of In Mz in [3] is a "spurious" result of the naive 
constituent quark model. If the conclusion of [6] were correct, it would call 
into question the validity of QCD studies in a variety of contexts, since it 
suggests that low energy strong dynamics may influence very short distance 
phenomena, a violation of the basic tenet of asymptotic freedom. 

On the other hand, with respect to large distances, the analyses in [5,6] 
correctly emphasize differences in the hadronic dynamics of longitudinal 
and transversal parts of quark triangles. These differences are not properly 
accounted for by approximation with constituent masses of light quarks. 
Moreover, constituent masses break the chiral structure of the fermion loops. 

The above challenge prompted us to reanalyze the triangle contributions 
to g — 2 in [7]. We believe to have located the source of the error in [6]; 
our arguments are summarized in this text. We also address some questions 
raised after [7] was published. Let us stress that the discrepancy in question 
is not relevant for an interpretation of the present g — 2 experiment, whose 
design accuracy is ±40 • 10 -11 . This is an academic dispute par excellence. 
However, its clarification is important, especially if our community is to 
determine the much more difficult hadronic three-loop light-by-light effects 
that limit the accuracy of the standard model prediction for g — 2. 
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2. Hadronic triangle diagrams 

Some details of the virtual fermion triangle are shown in Fig. 2. For the 

external magnetic field 




f 

Fig. 2. Virtual fermion loop with structure of couplings and external momentum. 

determination of the muon anomalous magnetic moment, we are interested 
in the Z* — ► 7* transition between virtual Z and 7 in the presence of the 
external magnetic field to first order in this field. Moreover, the magnetic 
field is constant so Z* and 7* carry the same momentum q. In this ap- 
proximation, but including all effects of strong interactions in the quark 
loop (such as gluon exchanges and non-perturbative effects), the amplitude 
of the Z* — 7* transition T fiu (q) can be parametrized using two Lorentz- 
invariant functions of the external momentum q, 

T M = ^2 [wT(q 2 )(-q 2 F^ + q^q a Fau-quq a F^) + w L (q 2 ) q v q° ' F aix ] . (3) 

Here, F^ v = (l/2)e flul sF lS denotes the dual of the external electromagnetic 
field tensor. The structure (3) is obtained in the following way. We construct 

which is a pseudo-tensor under Lorentz transformations using the vector 
q„ and the pseudo-tensor F^ u . From these objects we can construct three 
structures: F^ u , q^q a F au , and q u q a F afM . In addition, T^ v must be consistent 
with electromagnetic current conservation, q^T^ u = 0. This leaves us with 
two independent structures and we group them in such a way that one is 
transversal and the other longitudinal with respect to q v . 

The existence of the longitudinal part stems from the axial ABJ anomaly 
[9] which in fact exactly determines wi{q 2 ) in the chiral limit of massless 
fermions /, 

<f T ^ = ^L(q 2 )qVF^. (4) 
Comparing this with the axial anomaly we find 

w L {q 2 ) = ~ ■ 2Ip f Q) ee A . 2ip f Ql (5) 

where Q f and J| are the electric charge and the third component of weak 
isospin projection while Nj accounts for the number of colors in the case of 
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quarks and equals one for leptons. There are no corrections to this result in 
the limit m/ = 0, neither perturbative — due to Adler-Bardeen theorem [4], 
nor nonperturbative — due to 't Hooft consistency condition [10]. 

There is no such strong result for wt, but there is a very useful theo- 
rem proven by A. Vainshtein [8]. It states that for massless fermions, the 
following relation holds to all orders of perturbation theory: 

2w T (q 2 ) = w L (q 2 ) (for m f = 0). (6) 



There are nonperturbative corrections to this result which are given by pow- 
ers of Aq C y)/Q 2 but no perturbative ones. One way to prove Vainshtein's 
theorem is to consider the imaginary part of T„ v . The crucial point is that 
ImTp is symmetric under /i <-> v, q <-> —q. Indeed, since Im T^ u is given by 
convergent diagrams, we can freely use the anti-commutation of 75 to move 
it from the axial vertex 7^75 to the vector one, 7^. In the limit rrif = 0, 
this involves commuting 75 with an even number of 7 matrices, no matter 
how many gluon emissions occur on the way. 
Thus, we find 



Im 



w T (q 2 ) {-q 2 F^ v + q^q a F av - q u q a F a ^ + w L (q 2 ) q v q a F apu ^ 
^q^q a F au + q u q a F afM . (7) 



This is possible only if 

2ImwT = Imwi = const5(q 2 ). (8) 



Only with such a form for lmwT,L(q 2 ), Im?)^ has no terms antisymmetric 
in ji <-> v. The coefficient of the delta-function in Eq. (8) is fixed to be 
2-7T • 2ljNfQ 2 f by the exact form of wl in (5). Both wr(q 2 ) and u>L{q 2 ) 
are analytical functions decreasing at large q 2 . It means that they satisfy 
unsubtracted dispersion relations 

2 1 r Im w t ,l( s ) m 
7T Jo s - q z 

which when combined with Eq. (8) for imaginary parts implies relation (6) 
for the real parts as well. 

Note that the real part of T^ v contains an antisymmetric term F^ v . In 
terms of diagrammatic calculations it arises from counterterms fixed by the 
conservation of the vector current. When the Pauli-Vilars regularization is 
used these counterterms are given by heavy regulator loops. 
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wt = — —p^i — + non-perturbative corrections (10) 



To summarize, we now know that in the chiral limit rrif = 0, the longi- 
tudinal function wl is exactly given by Eq. (5) and the transversal function 
wt is known up to nonperturbative corrections, 

2ljN f Q 2 

~Q~ 2 

At this point we can apply our findings to an analysis of logarithmic con- 
tributions In Mz to g — 2. We will use the fact that Eqs. (5) and (10) rig- 
orously determine the large Q 2 asymptotics of wl,t', the non-perturbative 
corrections to wt, to be discussed below, are suppressed by extra powers of 
Aq C y)/Q 2 an d do not influence In Mz terms. Corrections due to deviations 
from the chiral limit are also suppressed by powers of rrif/Q. 

How do wl and wt enter the g — 2 contribution of the diagrams in Fig. 1? 
It was shown in [7] that for a determination of In Mz terms we can use the 
following simple representation: 




m| 



W L + ^2 _f n2 W T ■ (11) 



m 2 z + Q 

Two points can be made regarding this integral: 

• J°° dQ 2 WL diverges. The theory is inconsistent unless the anomaly 
cancellation condition is satisfied, 

J2 I f N fQ 2 f = ° (i 2 ) 

/ 

• /°° dQ 2 M z w T /(m 2 z + Q 2 ) ~ lnM z . With w T = w L /2 at large Q the 
anomaly cancellation condition (12) leads to the cancellation of lnM^ 
within a generation where rrif <C Mz for all fermions. 

In the analyses [5,6] the authors missed the leading l/Q 2 part in wt{u, d, s) 
but not in wl(u,cI,s) and wt,l(^, M> c )- Naturally, they arrived then at the 
expression (2) where In Mz terms do not cancel. Technically the mistake 
in [5, 6] arose from an incorrect construction of the Operator Product Ex- 
pansion (OPE) in the part related to wt- Referring to [7] for details note 
that the authors missed the leading operator which reflects interaction with 
the soft electromagnetic field (external magnetic field) at short distances. 
As we demonstrated above this leads to an inconsistency. 

On the other hand, the authors of [5, 6] were correct in their analysis 
of nonperturbative corrections, they were the first to point out that these 
corrections start with Aq CD /Q 6 in wt- So, below we will try to explain 
following [7] what can be done when we account for both the leading per- 
turbative term l/Q 2 in wt and subleading nonperturbative terms Aq CD / Q 6 
and higher. 
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3. Large Q 2 expansion, sum rules and models for Im wt,l 

For simplicity we limit ourselves in this section to the first generation, 
i.e. u and d quarks. As we mentioned above, the first nonperturbative 
corrections in the chiral limit m u ^ = are of order Aq CD /Q 6 in wt, as has 
been shown in detail in [6,7]. A particular example illustrating how the d = 6 
operators appear in the OPE is given in Fig. 3. The diagrams allow for a 




Fig. 3. Diagrams for four-fermion operators responsible for leading non-perturbative 
corrections to wt- 

perturbative calculation of the OPE coefficient while averaging of the four- 
fermion operator in the external magnetic field involves nonperturbative 
physics. 

Of course, the expansion of wt in powers of 1/Q 2 continues further 
and we are not able to find the sum. Still the analysis gives some exact 
relations: the first term is 1/Q 2 with the coefficient 1 (it is normalized to 

d^HN fQ 2 = 1), while the coefficient of 1/Q 4 is zero. These relations 
can be rewritten in the form of sum rules for Im.WT- Indeed, using the 
dispersion representation (9) at large Q 2 we find 

/•oo rco 

/ ds Imiuj'(s) = 7r , / ds s Im«ir(s) = . (13) 

JO JO 

The existence of such sum rules has been pointed out in [11], where their 
exactness was considered puzzling. It does not look puzzling to us, as there 
are many examples of such exact sum rules like Weinberg sum rules or sum 
rules for the DIS structure functions. Moreover, even if wt would not have 
the leading 1/Q 2 term, as assumed in [6], these sum rules would be no less 
exact (except the right hand side would be zero in the first one). 

It is instructive to compare these relations with similar sum rules for the 
longitudinal function wl where all powers of 1/Q 2 are known, 

roc /*oo 

/ dsImw L (s) = 2ir, / ds s k hnw L (s) = (k = 1, . . . , oo) . (14) 
JO JO 

This set of relations implies a unique solution Im^(s) = 2it5(s) showing 
that the massless pion is the only contributing intermediate state. 
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For the transversal function wt the intermediate hadronic states have 
to be 1 + mesons with isospin 1 and or 1 _ mesons with isospin 1. The 
lightest ones p, uj and a\ are massive even in the chiral limit. Representing 
ImwT as 

lmw T = ir^2gi5(s - m 2 ) (15) 

i 

we get from Eq. (13) 

i i 

The analysis of [6] assumed Y^iQi = 0- 

Until now we did not make any specific assumptions. Let us now assume 
a saturation by the lightest states and consider p and u to be degenerate. 
Then the solution of relations (16) gives 



m 



wt 



m 



ai 



— mt 



Ol 



) 2 + ™\ 



(17) 



This model for wt by construction satisfies exact relations (13) and was 
used in [7]. Of course, one can modify the model adding extra resonances 
but numerically is not much sensitive to these modifications. 

One more comment on the chiral symmetry breaking by quark masses. 
Accounting for nonvanishing m U) d leads to negligible effects for wt- The 
main (and important) effect for wl is a shift of the pole from to m\, i.e. 
to w L = 2/{Q 2 +ml). 

These refinements of wt,l in the long distance hadronic triangle effects 
for the u, d quark loops, together with similar ones for the s quark, were 
made in our study [7]. Numerically it changes the result in Eq. (1) to 
Aa^ MV {e,u,d;fi,c,s) = -6.7(1.0) • HT 11 . 



4. Summary 

Although the combined hadronic and leptonic triangle loop effects in 
electroweak corrections are not that significant for the muon g — 2, it is 
reassuring that a theory based on the Operator Product Expansion and its 
matching with hadronic phenomenology allows for quite accurate calcula- 
tions. We resolved a conceptual controversy existing in literature concerning 
cancellation of short distances between fermionic triangles. It seems to be 
possible now to develop an analogous approach to improve theoretical ac- 
curacy for the light-by-light part of hadronic effects. 



talk printed on February 7, 2008 



9 



Acknowledgments 

This research was supported in part by the Natural Sciences and Engineer- 
ing Research Council of Canada and by the DOE grants DE-AC03-76SF00515 
and DE-FG02-94ER408. 



REFERENCES 



[1] E. A. Kuraev, T. V. Kukhto and A. Schiller, Sov. J. Nucl. Phys. 51, 1031 
(1990); 

T. V. Kukhto, E. A. Kuraev, A. Schiller, and Z. K. Silagadze, Nucl. Phys. 
B371, 567 (1992). 

[2] T. Kinoshita and W. J. Marciano, in Quantum Electrodynamics, edited by 
T. Kinoshita (World Scientific, Singapore, 1990), pp. 419-478. 

[3] A. Czarnccki, B. Krausc, and W. Marciano, Phys. Rev. D52, 2619 (1995). 

[4] S. L. Adlcr and W. A. Bardeen, Phys. Rev. 182, 1517 (1969). 

[5] S. Peris, M. Perrottet, and E. de Rafael, Phys. Lett. B355, 523 (1995). 

[6] M. Knecht, S. Peris, M. Perrottet, and E. de Rafael, JHEP 11, 003 (2002). 

[7] A. Czarnecki, W. J. Marciano, and A. Vainshtein, Phys. Rev. D67, 073006 
(2003). 

[8] A. Vainshtein, Phys. Lett. B569, 187 (2003). 

[9] S. L. Adler, Phys. Rev. 177, 2426 (1969); J. S. Bell and R. Jackiw, Nuovo 

Cim. A 60, 47 (1969). 
[10] G. 't Hooft, in Recent Developments In Gauge Theories, Eds. G. 't Hooft et 

aL, (Plenum Press, New York, 1980). 
[11] E. de Rafael, talk given at the International Symposium on Lcpton Moments, 

Cape Cod, June 2003, http://g2pcl.bu.edu/~leptonmom/talks/deRafael.pdf, 

and private communication. 



